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I. INTRODUCTION 



Precision studies of top quark properties at the Fermilab Tevatron pp and CERN LHC pp 
colliders keep probing the Standard Model (SM) of electroweak (EW) and strong interactions 
at an increasing level of precision, and may provide a window to new physics. While the total 
ti cross section and tt invariant mass (M^f) distribution agree with SM predictions within their 
respective uncertainties, a measurement of the corrected (parton-level) forward-backward 
asymmetry in top-pair production, Apg, by the CDF [1-3] and DO [4, 5] collaborations at the 
Tevatron 

AJ?b(CDF[3]) = 20±7stat±2syst% 
A^b(D0[5]) = 19.6 ±6.5% (1) 

differs by about 2a [3, 5] from the SM QCD prediction. The difference between measure- 
ment and SM QCD prediction is even more pronounced, i. e. at the 3a level, in the region 
Mff > 450 GeV, where a measurement of Apg yields [2]: 

A^?b(CDF,M„- > 450 GeV) = 47.5 ± 1 1.4% (2) 

The forward-backward asymmetry in the ?f center-of-mass (CM) frame is defined as: 

.a ^ g/f (Ay > 0) - otii^y < 0) ^ 
™ arr{Ay>0) + arr{Ay<0) 

where Ay = yt —yt denotes the difference in rapidity of the top and anti-top quark, and SM 
predictions including higher-order QCD and EW corrections are provided in Refs. [6-13] 
and Refs. [12-14], respectively. The interpretation of the observed discrepancy requires a 
solid understanding of the theory predictions, i. e. control of the theoretical uncertainties. 
Recent updated calculations and studies of theoretical uncertainties at next-to-leading order 
(NLO) and next-to-next-to-leading-logarithmic order (NNLL) in QCD [11, 15] and at NLO 
EW-i-QCD [12, 13] find that the discrepancy is reduced compared to NLO QCD predictions, 
but still persists for the measurement at large Mft- For instance, including the EW contribu- 
tions to iff{a^) and ^{aaj) to A^p^g resuks in a combined NLO QCD-i-EW prediction of [12] 
(including the factorization/renormalization scale uncertainty and using MRST2004QED) 

A^B(NLOQCD-hEW) = 8.93+°;^^% 
A^B(NLOQCD-i-EW,Mff >450GeV) = 12.77+J;J^%, (4) 

and when including NNLL QCD contributions Apg is predicted as [11] (including the fac- 
torization/renormalization scale uncertainty and using MSTW2008): 

1 + 1. 04c 



A'p'b(NLO-i-NNLLQCD) = 7.24+^;^^% 

hl.7c 
-0.9' 



Aj[B(NLO-i-NNLLQCD,M„->450GeV) = 11.1+J;J%. (5) 
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However, since the first non- vanishing contribution to Apg is of NLO in QCD in the tt pro- 
duction cross section, a conclusive answer concerning the theoretical uncertainty will only 
be possible once a calculation of the complete NNLO QCD corrections to tt production 
becomes available. Nevertheless, it is interesting to study the possibility that the observed 
discrepancy could be interpreted as a signal of new physics. Possible SM extensions which 
may give rise to large contributions to Apg have been explored extensively in the literature 
and some recent examples can be found in Refs. [16-19]. In this paper, we consider the 
one-loop &{as) SUSY QCD and (ff{a) SUSY EW corrections to the strong partonic rf pro- 
duction process, qq tt, and study their impact on Apg at the Tevatron. We work in the 
complex Minimal Super symmetric SM (MSSM) [20, 21] and assume it to only be restricted 
by the condition of minimal flavor violation (MFV) [22, 23]. The other possible partonic 
process at leading-order (LO) QCD, gg — )■ tt, is symmetric in the production rates for top 
quarks in the forward and backward hemisphere, and thus only enters the total cross section 
in the denominator of Apg . After a consistent perturbative expansion of Apg in tts and a, the 
MSSM one-loop contribution to Apg calculated in this paper can then be written as follows 
(the dependence on and a is explicitly shown): 

A'^ = 0Cs^^ + cc^f-- (6) 

with Aa^QCD,SEW ^ S^SQCD,SEW^^y ^ 0) _ S^SQCD.SEW ^ ^^^^^ ^(0) ^^^^^^^ 

total tt production cross section at LO QCD and 5cr;^2C£»,5£H/ ^^^^^^ SUSY QCD and 
SUSY EW one-loop contributions, respectively. The SUSY QCD and SUSY EW one- 
loop corrections to (unpolarized) ti production in hadronic collisions have been studied 
in Refs. [24-31] and Refs. [31-35], respectively. They are known to only modestly im- 
pact the total tt production cross section and invariant tt mass distribution, and thus, at the 
presently available precision, do not spoil the good agreement between theory and exper- 
iment for those observables, at least for sparticle masses of ^(100) GeV (and larger) and 
rug > 230 GeV [30]. In this paper, we will derive general bounds on Apg, i. e. we work 
within the MFV-MSSM and scan over a large range of values for the relevant MSSM input 
parameters, without imposing additional constraints. In particular, these bounds do not rely, 
for instance, on specific SUSY breaking scenarios or an artificially reduced parameter space 
such as the constrained MSSM (CMSSM). As will be discussed as well, the impact of any 
specific assumption, e. g., mass limits derived from LHC squark and gluino searches, on the 
bounds on Apg can then be readily deduced from these general results. 
The paper is organized as follows. After presenting the calculation of the MSSM one-loop 
corrections that contribute to the forward-backward asymmetry in Section II, we derive an- 
alytic expressions for bounds on Apg induced by SUSY QCD one-loop corrections in Sec- 
tion II A and discuss the structure of the SUSY EW one-loop corrections to Apg in more 
detail in Section II B. In Section III, we present numerical results for the forward-backward 
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asymmetry. In Section III A, we determine general bounds on Apg based on the analytic 
expressions derived in Section II A. In Section IIIB, we present a comprehensive scan over 
the relevant complex MFV-MSSM parameter space. We conclude in Section IV and pro- 
vide explicit expressions for the relevant MSSM couplings and one-loop corrections in the 
appendix. 



II. MSSM ONE-LOOP CONTRIBUTIONS TO A 




Figure 1 : Momentum and mass assignments for the direct and crossed box diagrams. A possible color 
factor has been omitted. 

In the MSSM the additional one-loop contributions to the top forward-backward asymmetry, 
i. e. Apg of Eq. (6), originate from box diagrams involving squarks, gluinos and neutrali- 
nos inside the loops. The generic diagrams are shown in Fig. 1. The external momenta are 
denoted by pi , . . . , />4 and mi , . . . , m4 are the masses of the internal particles. The gluino- 
squark-quark and neutralino-squark-quark vertices are represented by r„ (n = 1,...,4), 
which we decompose, excluding a possible color factor, in the following way: 

r„ = g+p++g-p- with P^ = ^^ ■ (7) 

The full gluino-squark-quark vertex also includes color matrices FJ, = T"^ ■ r„ with T"^ = 

jl and the Gell-Mann matrices. The internal masses m„ and the generic couplings 
depend on the particles that are assigned to the internal lines. Explicit expressions for 
for gluino-squark-quark and neutralino-squark-quark couplings are provided in Section V A 
and Section VB, respectively. For each assignment of internal particles to the direct box 
(Fig. la), we get a corresponding contribution from the crossed box (Fig. lb) with identical 
values of the internal masses m„ and the coupling parameters g^. The vertices r„ (n = 3, 4) 
are related to the r„ by 

The amplitudes M^^^ and M^^^ corresponding to the diagrams in Fig. 1 can be generated 
using the Feynman rules for fermion number violating interactions of Ref. [36] that have been 
implemented in FeynArts [37-39]. They are then calculated with standard trace techniques 



4 



and a Passarino-Veltman reduction [40] of the loop integrals using FORM [41] . For the direct 
box diagram we write the partonic differential cross section in terms of the interference 
of M'^^^ with the LO matrix element as 



dcos G 



32;rf 

A 

32ns 



•A^, £ 2Re gt'g2'8fgt'D+-++ 

A = i 



(f, COS0) 

(f, COS0) 
(f, COS0) 
(f,COS0) 



(8) 



with Pt — y I — ^'^d denoting the top quark mass. The sum over the index A has to 
be interpreted as 



with ^± defined in Eq. (28) for SUSY QCD and in Eq. (31) for SUSY EW loop diagrams. 
The factors 1/4 and 1/9 in Eq. (8) are the spin and color average factors, respectively. The 
functions D^{s, cos 6) with 



a e {+-++,- + ++,- + +-,- + -+} 



(9) 



are given in Eq. (33). They depend on the partonic Mandelstam variables s={pi+ = 
{P3+P4)^ andby f = {pi —p^Y = {pi—PaY on the cosine of the scattering angle Q, i.e. the 
angle between the spatial components of p\ and p^,. They also depend on the internal masses 
rui, but these arguments are suppressed in Eq. (8) and the following discussion. Also, in 
Eq. (8) (and in the following discussion) it is implicitly understood that one has to sum over 
all possible assignments of particles to the internal lines. Note that the functions D'^{s, cos 0) 
do not depend on the generic couplings g^. The factor Na is the color factor and contains the 
factor 1/9 from the color average. The factor 1/4 from the spin average is absorbed into the 
functions D"^{s, cos d). 

The interference of the crossed box (Fig. lb) with ^^^^ can also be expressed in terms of the 
functions D"(5,cos 0). To see this, we note that the crossed box is obtained from the direct 
box by interchanging the momenta p^, and p^, reversing the fermion flow on the outgoing 
legs and replacing Yj, and r4 by r3 and respectively. The interchange of p^ and p^ 
is achieved by replacing cos0 by — cos0. When calculating the interference with the LO 
diagram, the appearance of fermion-number violating vertices must be handled correctly, 
e. g. following the rules in Ref. [36], and leads to an overall minus sign. The replacements 
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— >■ and r4 — ^ r4 are equivalent to — > ^J* (n = 3, 4). Thus, we obtain 
Jd(^) pt 1 1 



JCOS0 32;rf 4 9^^,^^ 

-g^^gt^g^^*g,^*D-+++{s,-cose) 
-gi^gt^g3^*gt^*D-++-{s -cose) 
-gi^gt^gf*g4^*D-+-+{s-cosG)]. (10) 

with of Eq. (33). Note that the interference of the crossed box with the LO matrix element 

may have a different color factor than the direct box. This factor is denoted by Nt,- 
Using Eqs. (8), (10) and do = dd'^"^ +dd^'^\ the corresponding partonic forward-backward 
asymmetry can be written in a compact form {d^j^ denotes the partonic LO total tt cross 
section, including both the qq and ^^-initiated tt production processes) 

ut 1 /"^ , .^^^^(^',cos0) /-o , Ja(f,cos0) 

^FB * = ^HT / ^^COS y ^ - / ^^COS : ' (11) 

{t^^Uq ^ ' dco^e .1-1 ^ ' dcosO ^ ' 

1 ,/ r^a(5,cos 0) ^^d(f, — COS0), 
= —7^ / a(cos0 — ^ — ^ ^- — -\ 

tt 

= /'j(cos0)£7?e{G^-A''(f,cos0)} 



with the index a defined in Eq. (9), 



A^(5,cos0) = -y^■A^g■[2D°(5,cos0)-2D^(5,-cos0)] (12) 



and 



= Lx=+,-gfg2^ (H^gt^Na + g^^*g,^*Ni,) , 
+++ = U=+.-gT^g^^ {gt^gt^Na + gJ^*g,^*Nh) , 
G,++- = U=+,-gl^gt^ igt^g4^Na + gi^*gf*Nt) , 

=U=+^-gi^gtHg3^gt^Na + g+^*g,^*Nh). (13) 

Note that the coupling products may, in general, depend on the initial-state quark flavor 
q, as indicated by the subscript. The color factors A^^ ^ have been absorbed into the definition 
of and all the other factors into the definition of cos0). An additional coupling 
factor Ng has been absorbed in A" for convenience so that the coupling products G^, defined 
in terms of |^ of Eqs. (29) and (32), are of ^(1). 
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The hadronic forward-backward asymmetry is then defined by folding with the parton distri- 
bution functions (PDF) and by dividing by the total hadronic LO cross section , so that 
A^g of Eq. (6) at the Tevatron reads 



A'' 



FB 



1 



^ / dxi / dx2 / d{cos6 



a 

+ /^/;,(^i)/,/p(^2) Y.^e{G\- - COS 0) } 0(f - 4m2) 

a 



(14) 



where we defined for each initial- state quark flavor q 



dxi dx2 d{cose)[fq/p{xi)fq/p{x2)-fq/p{xi)fq/p{x2)] 



a(«) ^0 JO JO 



X G{s-4mj)-A'\s,cosd) 



and 



Re{ G, 



+ Gq^^-Aq 



(15) 



+ Gq+-+Aq 



} 



(16) 



Here -y/S = 1.96 TeV is the Tevatron hadronic CM energy. The function fq/p{fq/p) is the 
PDF of the quark(anti-quark) flavor q{q) inside the proton and the functions fq/p and fq/p 
are the corresponding PDFs for the anti-proton. In Eq. (15) we made use of the fact that 
at a pp collider the quark( anti-quark) distribution inside the proton coincides with the anti- 
quark(quark) distribution inside the anti-proton. Moreover, when the incoming momenta are 
interchanged (or xi ^ X2), it corresponds to replacing cos with — cos 6. The difference in 
the PDFs occurs because A'' of Eq. (12) is anti- symmetric in cos0. 

In the presence of CP violating phases, the coupling factors Gq of Eq. (13) may have imagi- 
nary parts. In this case, the imaginary parts of the may also contribute to Apg. 
With the notations above, we have decomposed the MSSM one-loop contribution to the 
forward-backward asymmetry, Apg of Eq. (6), into loop functions A^ and coupling products 
Gq. The loop functions only depend on the internal masses in the box diagrams of Fig. 1 
while the dependence on the coupling constants is contained in the coupling factors G°. This 
separation will prove very useful when we attempt to give bounds on Apg in the MSSM that 
do not rely on specific SUSY breaking scenarios or an artificially reduced parameter space. 
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A. SUSY QCD One-loop Contributions 



When only SUSY QCD one-loop contributions to Apg are considered, the Majorana fermions 
in the box diagrams of Fig. 1 are all gluinos. The relevant diagrams are shown in Fig. 2. In the 




Figure 2: Feynman diagrams of the SUSY QCD one-loop contribution to the forward-backward 
asymmetry Aj^g. 

case of SUSY QCD, Eq. (16) can be further simplified due to relations among the coupling 
products and the loop functions A^. The color factors in Eq. (13) and the coupling factor 
in Eq. (12) are now 

7 2 



where gs is the strong coupling constant. The internal masses are 

m2 = m4 = rUg , mi = triq. , mj, = mf. , 

where / and j are the sfermion indices of the squark and the stop, respectively. To make the 
dependence on the sfermion indices explicit we use the notation 

Gq Gq ij , A'^ ij. 

Ultimately, we must sum over the sfermion indices i and j. However, since both the internal 
masses and the couplings may depend on the indices i and j, we postpone that sum until the 
end of the discussion. Due to m2 = m4, we have the relation 

A+-++ - 4- + + + 

as can be seen from the expressions for D'^ in Section V C. 

Using the unitarity of the squark mixing matrices of Eq. (30), we can simplify sums and 
differences of the coupling products: 

+ = -2{Na+Nb)Re{U]^^u'*_e'^) ^ {Na+Nh)Gf^ , 

Thus, we write for A„ of Eq. (16) (again we added the sfermion indices ij) 
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with 

4)j = iNa+N,)^{A^±^-+A;^+r+) (18) 

and 

Cf^ = -2Re{ul+U^*_e"f) 

Using the unitarity of the stop mixing matrix U\ we find 

G<" = -C«" , GS = -Gi^ . C^X = -cgr (20) 

Note that G^p and G^^-j are always real, and one therefore only needs to consider the real 
part of the functions -j. Furthermore, the coupling product G^ does not depend on the 
initial-state quark flavor q. If one assumes that the mixing matrix of the light flavor squarks 
q is diagonal, as in almost all of the considered parameter space, G„ / ■ reduces to 

G^j = {-iy-{\ui^-\u[_\')^Gf (21) 

and does not depend on the initial-state quark flavor q either. We can thus perform the sum 
over /, j and q and obtain 

<sqcd = gS^^a(i)+gS?a(2)+a(3) (22) 

with 

q=u,d,s,c 
q=u,d,s,c 

^"'= E ''-K'!.+4'i2+4i+4ii. (23) 

q=u,d,s,c 

Note that A^^^ vanishes for degenerate stop masses and A^^^ vanishes if the stop or squark 
masses are degenerate. Upper and lower limits on the SUSY QCD one-loop contribution to 
Apg can now be obtained for given values of the squark masses by exploiting the fact that the 
limits of g[^'^ and G^"* are 

-1 < g\^^ < 1 , -1 < Gjf < 1 

and thus 

^(3)_|^(l)|_|^(2)| <aJ[b,sqcd<^^^^ + I^^^^I + I^^^^I ■ (24) 
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B. SUSY EW One-loop Contributions 



The SUSY EW one-loop contributions to the forward-backward asymmetry Apg of Eq. (14) 
consists of four diagrams, two direct boxes and two crossed boxes, shown in Fig. 3. We again 



IP 



IP 




U S t U Xm t 

(a) {b) 

Figure 3: Feynman diagrams of the SUSY EW one-loop contribution to the forward-backward asym- 
metry Apg. 

start with the generic expression for Apg in terms of loop functions and coupling products as 
given in Eq. (16). The color factors in Eq. (13) and the coupling factor in Eq. (12) are now 



4„2 



where e is the electromagnetic coupling. The internal masses are 

mi = m^;, m3 =mf., 
Fig.3(a,c): m2 = m^^, = trig or Fig. 3(Z7,(i) : m2 — mg, — m^,^ 

where / and j are the sfermion indices of the squark and the stop, respectively, and k is the 
neutralino index. To make the dependence on the sfermion and neutralino indices explicit, 
we use the notation 



Ga . /-•a 
q ^q.kij 



A" A" 



Ultimately, we must sum over the indices i, j and k. This are 16 combinations for each 
diagram. By examining the amplitudes for the individual diagrams one notices that the 
diagrams (a) and (c) in Fig. 3 are the complex conjugates of the diagrams (b) and (d), respec- 
tively. Thus, the are real and only the real parts of the coupling functions G° of Eq. (13) 
contribute. The contribution of the SUSY EW one-loop corrections to Apg of Eq. (16) then 
reads 



^FB,SEW 



4 2 



q k=U,j=l 

+ Re{G 



q,kij 



-+ 



(25) 
(26) 
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In case of SUSY EW one-loop corrections, it is quite complicated to analytically derive 
bounds on the coupling factors , ■ because of the complicated structure of the neutralino- 
squark-quark couplings (see Eq. (32)). It is even harder then to find reasonable bounds on 
^'fb sew SUSY QCD case. The box diagrams have 16 different squark und 

neutralino mass combinations where always some cancellation occurs due to the unitarity of 
the mixing matrices. Therefore, we performed a MSSM parameter scan to extract bounds on 
A'lf, c-Fw as described in Section IIIB. 



III. NUMERICAL RESULTS 



For the numerical evaluation of the forward-backward asymmetry Apg of Eq. (14), we use 
the LO PDF set CTEQ6L1 [42] with the renormalization (iir) and factorization (jip) scales 
chosen to be equal to the top quark mass, iir = iif = mt. The SM input parameters are = 
173.2 GeV, a = 1/137.036, Mw = 80.36 GeV, Mz = 91.187 GeV, and cosdw = Mw/Mz. 
We assume one-loop running of the strong coupling constant with UsiMz) = 0.130, so that 
OCs{mt) = 0. 1 18, which is consistent with our choice of PDFs. To evaluate the coefficients of 
the tensor integrals Dj ij of Eq. (33), the LoopTools library [43] has been used. 
The study of the dependence of Apg on the MSSM input parameters is simplified by the 
fact that it is not sensitive to all parameters of the complex MFV-MSSM. First of all, the 
forward-backward asymmetry can only come from diagrams with up or down quarks in the 
initial state. For strange, charm and bottom quarks, the PDFs are the same as those of the 
corresponding anti-quarks so that the difference of PDFs in Eq. (15) is zero. In the MFV- 
MSSM Apg is therefore insensitive to parameters that only affect the masses and couplings 
of strange, charm and bottom-squarks. Furthermore, the trilinear couplings of up and down- 
squarks only enter through the squark mass matrices, where they are suppressed by the small 
up and down-quark Yukawa couplings. Thus, only the following set of MSSM input param- 
eters are relevant to our study: 

• tan/3 

• M 

• the pseudoscalar Higgs mass niA 

• the gaugino masses Mi, M2 and M3 

• the top- squark trilinear coupling At 

• the soft masses rriq^^ and rriq^^ of the left-handed first and third generation squarks 

• the soft masses m^^, mj and m^^ of the right-handed up, down and top squarks 
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Of these parameters, /i, Mi, M2, M3 and At can be complex, but one of these phases can be 
rotated away. We rotate the phase of M2 away and study the dependence on the remaining 
complex phases and the absolute values of the above MSSM parameters independently. 
These input parameters are constrained by direct SUSY searches at LEP, Tevatron and LHC, 
and indirectly by low-energy precision observables. A review of results from the search for 
signals of low-energy SUSY at LEP and the Tevatron as well as in precision observables 
can be found, e. g., in Ref. [44]. Most recently stringent exclusion limits on squark and 
gluino masses within the CMSSM and Simplified Models have been obtained at the LHC by 
the CMS [45-47] and ATLAS [48-51] collaborations. In Section IIIB we provide general 
upper and lower bounds on Apg within the complex MFV-MSSM by performing a compre- 
hensive scan over a wide range of values for these input parameters without any additional 
assumptions. In particular, these general bounds do not rely, for instance, on specific SUSY 
breaking scenarios or an artificially reduced parameter space such as the CMSSM. The re- 
sults are presented in such a way that the effect of a change in the sparticle mass limits on 
the upper and lower bounds on Apg can be estimated. In Section III B this is done for the 
current LHC squark and gluino mass limits. 

In the next sections we first discuss the main characteristics of the MSSM one-loop correc- 
tions to Apg and then present results of a comprehensive parameter scan. 



A. SUSY QCD and SUSY EW Contributions to the Loop Functions 

The main characteristics of the MSSM one-loop corrections to, and of bounds on, Apg can 
be determined by studying the SUSY QCD loop functions a'^^I of Eq. (18) and SUSY EW 
loop functions -^^ of Eq. (26). Here we will only present results for the contribution of the 
MM-initiated tt production process, since the dd production channel is much smaller due to 
the smaller PDF (about a factor of eight smaller as discussed in Section IIIB). 
In case of SUSY QCD contributions we showed in Section II A that bounds on Apg sqcd 

(k) 

Eq. (18) can be derived in terms of the normalized hadronic loop functions A^ jj.k = 1,2, 3, 
where q denotes the initial-state quark flavor, / the squark index / = 1,2 with flavor q, and 
j = 1,2 refers to the top squark index. These loop functions only depend on three general 
mass parameters, mi = rtiq., m2 = nig, and mj, = mf., as indicated in Figs. 1,2. For the 
special case of initial-state up-type quarks and mass degeneracy, m\ = m2 = m^, = M, we 
show numerical results for the functions A^^ ■ ■ {k = 1,2,3) in Figs. 4(a) and 5. As can 

' (2) 

be seen in Fig. 4(a), the largest single contribution to the loop functions comes from A^^ -j. 
This feature persists when the masses are varied independently. A^^ '■■ reaches up to + 0.9% 
for M ^ 100 GeV and - 0.5% for M ^ 200 GeV. The other two loop functions can reach 
roughly -|- 0.4% for M ^ 100 GeV and only tiny negative values. For masses M > 400 GeV 
all contributions are very small. The first peak in Figs. 4(a) and (b) at M = 86.6 GeV is due to 
a normal threshold [52] when the condition p\ = {1712 + m^)^ or = (ms + m4)^ is fulfilled. 
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Figure 4: (a) SUSY QCD normalized hadronic loop functions A\ [. with k = 1,2,3 as defined in 
Eq. (18) and (b) bounds Amm,max on A^g sqcd ^1- (^^) ^^^^ initial-state up-quarks. Shown is the 
dependence on a common mass M. The bounds are obtained assuming large up-squark and top- 
squark mass splittings. No kinematic cuts have been applied. 



which here is the case when M = 7%/2. A resonance in the partonic tt cross section occurs 
when the gluino pair in the box diagrams can become on-shell, thus nig > nit. This resonance 

(k) 

manifests itself as dips in the hadronic functions AJ^ '-j at M ^ 200 GeV. When the resonance 
occurs inside the numerically important tt invariant mass region of = 350 to 500 GeV, it 
leads to a correspondingly large negative asymmetry. For larger SUSY particle masses M, the 
resonance region is shifted to larger values of and therefore outside the dominant region 
of the tt cross section. From these results an estimate of the bounds A,nin,max on A^g sqcd 
can be obtained using Eqs. (23) and (24). As noted earlier, A^^^ vanishes for degenerate stop 
masses and A^^^ vanishes if the stop or squark masses are degenerate. In general, we found 
that the SUSY QCD one-loop corrections to the forward-backward asymmetry increase if the 
up-squark mass splitting or, in particular, the top-squark mass splitting are increased. Thus, 
the largest asymmetry is obtained if and t2 are decoupled, which results in vanishing 

(k) (k) (k) 

functions 21 ' ^11 ■ '•'^^^ scenario the bounds of Eq. (24) read: 



Amin 
AvYinx 



.(3) 
^(,11 



\A^'^ I 



lA^') I 



pp') , I + IA^^) I 

^uA\ + l^«,lll + I^K,lll • 



(27) 



min,max 



on M = m? 



nig. The 



Figure 4(b) shows the dependence of these bounds A 
total asymmetry can reach values from — 0.9% up to + 1.6% for this configuration when no 
kinematic cuts have been applied. 

(k) 

Figures. 5(a) and (b) show the loop functions A^^ ■■ in dependence of a cut on Mff, M™'", and 
on |A_y|, Aymin, respectively. FromM^"' 



'J 

350 GeV to M'^f" 



600 GeV, the functions A 
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Figure 5: SUSY QCD normalized hadronic loop functions jj witii ^ = 1,2,3 for initial-state up- 
quarks as defined in Eq. (18), assuming all particles in the loop have a common mass M = 100 GeV. 
Shown is the dependence on the cut on (a) the invariant tt mass M,f, Mfj > M""", and on (b) the 
rapidity difference \Ay\ = \yt—yt\, \^y\ > ^ymm- 



andAf,'',-; roughly double in size while. a[}}- increases by about 50%. If one increases M"/'" 
from 350 GeV to 450 GeV, the loop functions increase by a factor of 1.35 for a|/-^, 1.43 for 
A^^P-j and 1 .52 for A^^'^^-^. Thus, the bounds on the forward-backward asymmetry for the SUSY 
QCD contributions will roughly change by a factor of 1.4 if one applies a M'™" = 450 GeV 
cut compared to the case without such a cut. Applying a cut on \Ay\, as shown in Fig. 5(b), 
also increases the loop functions A^^j^^'^'^'*. For instance, when changing Aymin from zero to 
one. All I J increase by a factor of 2.00 and A^^ Ij by a factor of 2.25. 

We now investigate the loop functions of the SUSY EW one-loop contributions to the 
forward-backward asymmetry as described by A^g^^yy of Eq. (26). In this case, one 
has four different masses in the loop and four different loop functions A"^ -j^. with a E 

{H — h ++, — h H — , — I h}. If one plots these loop functions in dependence 

of a common SUSY mass M, one can produce a similar plot as in Fig. 4, just rescaled by the 
smaller coupling factor and different color factors. To illustrate the effect of the neutralino 
mass we therefore show in Fig. 6 the loop functions for the SUSY EW one-loop contribu- 
tions as defined in Eq. (25) in dependence of the neutralino mass, when assuming a common 
mass for the other sparticles in the loop: M = mg = mj. = ma^ = 100 GeV. The loop functions 
"^utjk^^ and A^^^^ give the largest contributions. They can amount to about ±0.04% and 
only slowly decrease with increasing neutralino mass, e. g., they are still about ±0.02% for 
relatively large neutralino masses m^o ^ 700 GeV. 

As stated in Section n B, it is difficult to obtain an analytic expression for bounds on A^^ ^^^^ 
like it was done in case of SUSY QCD-induced asymmetries. We refer therefore to the next 
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Figure 6: SUSY EW normalized hadronic loop functions -y^ with a G {H — — h 

H — , — I h} for initial-state up-quarks as defined in Eq. (25) in the case where the up-squark, top 

squark and the gluino have a common mass M = mg = mj. = = 100 GeV. Shown is the dependence 
on the neutralino mass niyO. No kinematic cuts have been applied. 

section where we perform a comprehensive scan over the relevant MFV-MSSM parameter 
space to determine the largest possible value of A^g induced by both SUSY QCD and SUSY 
EW one-loop corrections. 

B. Bounds on A^g from a MSSM parameter scan 

In phenomenological studies of the MSSM the large number of parameters is a common 
problem. Even if all parameters of the soft MSSM Lagrangian are assumed to be real and 
all flavor structures are assumed to be proportional to the SM Yukawa matrices (minimal 
flavor violation) we are left with 30 independent parameters [53]. Numerical discussions 
of observables within the MSSM are therefore often limited to constrained scenarios, where 
certain assumptions about the SUSY breaking mechanism are imposed, or even to individual 
benchmark points which are deemed "representative" in some sense. Here we describe in 
detail how we performed a comprehensive scan over the relevant parameter space of the 
complex MFV-MSSM to determine bounds on Apg. 

In the case of Apg, we are actually able to scan the full relevant parameter space of the 
complex MFV-MSSM. This is made possible by the separation of Apg into loop functions 
A^ and products of coupling parameters G^, which is described in Section II. Calculating the 
loop functions of Eq. (15) requires a numerical phase space integration and is therefore rather 
time-consuming. However, the A^ are relatively smooth functions of only four variables 
(namely, the internal masses mi, . . ..m^). The loop functions may therefore be calculated 



mg =m- =mg = 100 GeV \ 




— ■ — + - + + 

+ + + 

- + + - 

+ - + 
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on a four-dimensional "mass-grid" and linear interpolation can be used to obtain for 
other mass values within the grid range. The remaining computational cost of calculating 
masses and coupling constants for a given set of MSSM parameters and interpolating the 
loop functions is extremely small, so that sampling up to one billion (!) MSSM parameter 
points is perfectly doable on a single core computer. 

Let us briefly discuss the details of the grid interpolation and the parameter scan. Grid data 
for the loop functions were generated for masses between and 2 TeV. For masses below 
500 GeV the grid spacing was 20 GeV. If one of the masses exceeds 500 GeV the grid spacing 
was increased to 50 GeV in that direction. If a mass exceeds 1 TeV the grid spacing was 
increased again to 100 GeV. At each grid point the integrals in Eq. (15) were calculated with 
the VEGAS algorithm. Specifically, we use the OmniComp-Dvegas package [54], which 
facilitates parallelized adaptive Monte Carlo integration and was developed in the context of 
[55, 56]. A lower cut of 450 GeV on the tt invariant mass was applied throughout. Thus, all 
the results in this section are for the 'large nitf bin. Separate integrations of were done 
for each value of the superscript a (three for the SUSY QCD contributions and four for the 
SUSY EW contributions), but the values for different quark flavors q were determined with 
the same simulation. Adaptation was driven by the M-quark flavor, which always produces 
the largest value. The relative accuracy of the numerical integration was required to be below 
1%. 

The multivariate linear interpolation was done by successively using one-dimensional linear 
interpolation in each of the variables. If, during the parameter scan, a certain mass exceeds 
2 TeV the loop functions where it enters are assumed to be zero. In doing this, we neglect 
loop functions of the order of one permille. The discrepancies between the exact and the 
interpolated values of the loop functions are of the same order. 

As discussed earlier, the parameter scan is simplified by the fact that Apg is not sensitive to 
all parameters of the complex MFV-MSSM. Our choice of relevant MSSM input parameters 
is listed in the beginning of Section HI. The value of tan /3 was varied between 1 and 50 
and all parameters with mass dimension one were varied between and 3 TeV. The complex 
phases of jU, M\, Mt, and At were varied between and In. 

For the actual scan we used an adaptive method along the lines of Ref. [57]. The basic 
idea is the following: Instead of sampling all parameters with a uniform (or otherwise fixed) 
random distribution, one defines an importance function which quantifies the importance of 
a given set of parameters. Since we are interested in MSSM parameter points with large 
effects in Apg we used |Apg | as importance function. We then used VEGAS to compute 
the integral of the importance function over all the scan parameters. Adaptation guarantees 
that the "important" regions of the parameter space are sampled with a higher density. The 
OmniComp-Dvegas package [54] was used for the parameter scan too. A total of 8 ■ 10^ 
parameter points were sampled. Adaptation was done with 22 iterations. 
With the data from our scan, we can now show results for upper and lower bounds on Apg 



16 



as a function of any relevant MSSM input parameter. To do this, we simply bin the sample 
points with respect to that input parameter and determine the maximal and minimal value of 
A^g in each bin. 

Fig. 7 shows the upper and lower bounds on Apg as a function of the gluino mass. Shown 
separately are the contributions from the uU and J J-initiated tt production channels as well 
as from SUSY QCD and SUSY EW one-loop corrections, assuming M„- > 450 GeV. 
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Figure 7: The gluino mass dependence of the lower and upper bound on A'^^, |Apg | < A, due to SUSY 
QCD and SUSY EW one-loop contributions, with Mfj > 450 GeV. Separately shown ai^e the bounds 
when only including the uil and JJ-initiated production channels. 

The bounds shown in Fig. 7 have been obtained by assuming a lower limit on the top squark 
masses of m^j , > 100 GeV and a mass limit for the other squark flavor masses of m^j , > 
300 GeV. For the neutralinos no mass limit was imposed. As discussed earlier, the uu- 
channel generates the largest contribution to Apg due to the large PDF and is roughly a 
factor of eight larger than contribution from the JJ-channel. The SUSY QCD bounds peak 
around lUg ^ 270 GeV, when the contribution of the resonant gluino pair in the box diagrams 
coincides with the maximum of the LO Mft distribution. The absolute values of the upper 
and lower bounds are practically identical because for the applied squark limits the SUSY 
QCD loop function A^^^- is dominant (see Fig. 4), so that the SUSY QCD bound of Eq. (24) 
is well approximated by — |a(^)| < A^g ^q^^ < \A^^^. 

As shown in Fig. 7, the bounds on the SUSY EW one-loop corrections to Apg are much 
smaller than the SUSY QCD one-loop corrections, since they are suppressed by the smaller 
electroweak coupling. Furthermore, the bounds are not as symmetric as in the SUSY QCD 
case. In general we found the absolute value of the lower bound to be larger than the upper 
bound. The SUSY EW one-loop corrections to the dd-chmnel are not shown in Fig. 7 
as they are basically zero, being suppressed by both the ^^-quark PDF and the electroweak 
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coupling. Again, larger SUSY EW loop-induced asymmetries can be obtained when relaxing 
the constraints on the MSSM parameters. For instance, for sparticle masses below 50 GeV 
one can obtain A^^ = —0.4%. 

Combining the SUSY QCD and SUSY EW one-loop contributions as well as taking into 
account all ^^-initiated tt production channels, the lower and upper bounds on Apg in 
the complex MFV-MSSM, A,„,„ and A,nax, are shown in Figs. 8 and 9. In Fig. 8 we 




Figure 8: Bounds on the MSSM one-loop contributions to Apg with Mfj > 450 GeV: (a) Lower and (b) 
upper bounds on A^g (in percent) as functions of the lighter top squark mass and m„„„, which is 
the lightest other SUSY particle in the loops except for mf^_ , i. e. m„,„, = Min{mg , m^j , , ^ , niuj ^,m^^^}. 

show the dependence of the bounds Amin^max on the lighter stop quark mass mj^ and on 
nimin, which is the lightest other SUSY particle in the loop except for mf^, i. e. m,nin = 
Min{mg,mj^j2 3 4,mui2,m;jj2}- small masses of irij^ < 200 GeV and m„j,„ < 200 GeV, 
the upper bound on Apg is with up to +3% somewhat larger than the absolute value of the 
lower bound with —2%. This is because the SUSY QCD bounds, Eq. (24), are in general 
not symmetric and in this region the contribution from A^ / is non-negligible and positive as 
shown in Fig. 4. 

In general, one observes that smaller up-squark masses mu^ and top squark masses nif^ lead 
to larger MSSM one-loop contributions to the forward-backward asymmetry. The same 
is valid for the gluino mass when irig > 250 GeV. For nig < 250 GeV the asymmetry is 
dominated by the contribution of the resonant gluino pair. Since the SUSY EW one-loop 
corrections are sub-leading, the asymmetry has a very small dependence on the neutralino 
mass. Furthermore, the asymmetry is larger for a larger mass splitting of mfj and m^,. This is 
illustrated for the lower bound on Apg in Fig. 9(a) and the upper bound on Apg in Fig. 9(b), 
where A,nin,max are shown in dependence of mj^ — m-f^ and on nimin, which here is the lightest 
SUSY particle in the loops, i. e. m„„„ = Mm{mi^^,mg,m.j^^^^^,m^^^,m~^^^. 
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Figure 9: Bounds on the MSSM one-loop contributions to Apg with M,,- > 450 GeV: (a) Lower and 
(b) upper bounds on Ajjg (in percent) as functions of the stop mass splitting nif^ — and tUmin, which 
is the lightest SUSY particle in the loops, i. e. m^m = Min{mtj2,mg,m^j 234)^101 21 rnSi 

As can be seen, when the stop masses are degenerate, then even for tiny m,„(„ one finds 
—0.5% < Apg < + 1 %. For a stop mass splitting of mf^ — = 500 GeV and nimm > 300 GeV 
one can obtain bounds up to \Amin^max\ = 0.5%. Larger values for the asymmetry are possible 
when the stop mass splitting is even larger and m„„„ < 300 GeV A similar behavior, but not 
as pronounced, can be observed for the up-squark mass splitting ma^ — nia^ . For the bounds 
shown in Fig. 9(a) and Fig. 9(b), U2 is basically decoupled and ma^ ^ w„„„. In the case 
when the two up-squark mass eigenstates are degenerate and ma ~ m„„„ > 100 GeV, one can 
roughly divide the given limits of Fig. 9 by a factor of two to obtain the limits for the case 
of degenerate up-squark masses mu^ = ma^. For ma ~ m,„,„ < 100 GeV the limits decrease 
only mildly. 

As discussed earlier, SUSY searches at the LHC, interpreted in the CMSSM and Simplified 
Models, impose stringent limits on the masses of squarks and gluinos. An updated summary 
of the limits obtained by CMS and ATLAS can be found in Refs. [58, 59]. In the following 
we illustrate the impact of some of these limits on the bounds on Apg. For instance, ATLAS 
SUSY searches in events with jets and missing transverse energy interpreted in the CMSSM 
find that gluinos and squarks lighter than about 950 GeV [48, 59] are ruled out for nig = niq, 
while for niq ^ nig the gluino can be somewhat lighter, i. e. trig > 680 GeV [59, 60] for 
niq = Inig ( all limits are at 95% C. L.). Note that these squark mass limits do not apply for 
the third generation squarks due to slightly different production mechanisms for stops and 
sbottoms and the different decay pattern of top squarks. As can be seen in Fig. 8, in this 
scenario with ntg = m,„in = 680 GeV the upper bound on Apg can be at most A,nax ~ 0. 15% 
for mr, < 200 GeV. 
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Note that these mass limits depend on the assumptions that go into the extraction of the limits. 
For instance, when rtig ^ niq the SUSY signal is dominated by squark pair production in the 
gg — 7- qq channel and thus is proportional to the number of degenerate squarks. The signal 
cross section drops significantly if all but one squark are decoupled, which weakens the mass 
limits obtained under assumptions of degenerate squarks. The sensitivity of these SUSY 
searches also drops when the neutralino mass m^^o is increased. For example, as shown in 
Ref. [58] the CMS limit on the gluino mass obtained in a Simplified Model decreases from 
about 850 GeV to 400 GeV, if the neutralino mass is increased from GeV to nig — 200 GeV. 
In this scenario, the upper bound on Apg increases from about 0.1% (m„j;>; = 850 GeV) to 
A,„ax ~ 0.5% (mmin = 400 GeV) for mf, < 200 GeV, as shown in Fig. 8. 
Moreover, the extraction of mass limits on squarks and gluinos gets even more involved if 
the (light) squarks do not decay directly into the lightest neutralino. And finally, one might 
consider MSSM scenarios where the neutralino is not the LSP. The LHC SUSY searches 
of Refs. [45, 48, 61] do not strictly require a neutralino LSR They remain valid as long as 
squarks and gluinos decay dominantly into missing energy and one or two jets, respectively. 
A different type of analysis is required if colored SUSY particles are stable or quasi-stable. 
This can, for example, happen in gauge mediated SUSY breaking scenarios with a gravitino 
LSP and a squark next-to-LSP (NLSP). Stable or quasi-stable SUSY particles would form 
so-called i?-hadrons with specific detector signals. A search for i?-hadrons was presented in 
Ref. [62] and mass limits of 300 GeV and 600 GeV are given for (quasi-)stable squarks and 
gluinos, respectively. Again with smaller squark/gluino mass, the bounds on Apg increase 
in this scenario from about 0.25% (m„„„ = 600 GeV) to almost \A,nin,inax\ ~ 1% (/«min = 
300 GeV) for m^- < 100 GeV, as shown in Fig. 8. 

As illustrated on these examples, since our results for the bounds on the forward-backward 
asymmetry have been obtained with the only assumption that we work within the MFV- 
MSSM, the impact of a number of specific scenarios on the upper and lower bounds on Apg 
can be estimated. 

rv. CONCLUSIONS 

The study of asymmetries in tt production, such as the forward-backward charge asymme- 
try, parity violating asymmetries in polarized tt production, and spin correlations between 
the t and t, may provide a window to non-SM physics complementary to direct searches for 
non-SM particles. A recent measurement of the corrected (parton-level) forward-backward 
charge asymmetry in tt production at the Tevatron deviates from the SM prediction by about 
3 a in the region > 450 GeV. Provided the SM prediction is under good theoretical con- 
trol, this discrepancy may leave room for an interpretation as a signal of non-SM physics, 
and a number of non-SM scenarios have been proposed in the literature. In this paper we cal- 
culated the SUSY EW and SUSY QCD one-loop corrections to A^g within the MFV-MSSM 
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and derived general lower and upper bounds Am,„ , A,nax on Apg at the Tevatron by performing 
a comprehensive scan over the relevant MSSM parameter space. Since the MSSM one-loop 
corrections to Apg are dominated by SUSY QCD one-loop corrections, these bounds strongly 
depend on the squark and gluino masses and are much less affected by the neutralino mass. 
As expected, these loop effects are more pronounced for smaller sparticle masses in the 
loops and larger stop/squark mass splittings. For very small masses, mf^ < 200 GeV and 
>n,nin = Min{mg,mu, 2,m^^ ^} < 200 GeV, we find —2% < A^g < +3%, which is comparable 
in size to the SM EW corrections to Apg. This is the maximum possible SUSY loop-induced 
asymmetry that can be obtained within the MFV-MSSM at the Tevatron, with M„- > 450 GeV. 
It is important to emphasize that apart from working within the MFV-MSSM, no additional 
assumptions or constraints have been imposed in our derivation of the bounds on Apg and 
that from the general bounds presented in this paper an estimate can be obtained of how they 
change when assuming a specific SUSY scenario. For example, to illustrate the impact of 
squark and gluino mass limits obtained within the CMSSM at the LHC, we obtain from the 
results presented in Fig. 8 that for gluino and squark masses in the range 850 GeV- 1000 GeV, 
^max ~ 0. 1 % for a light top squark of m^^ = 200 GeV and /n,,,,-,, = 850 GeV, and A,,,^;^ ^ 0.05% 
for mf J = 600 GeV and m,„,„ = 1 TeV 
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V. APPENDIX 

A. Gluino-squark-quark couplings 

The gluino-squark-quark couplings as defined in Eq. (7) are [39] 

= iVlgsgt with (28) 
g"^ = Xuy^^l^ , |i = At/jy^^/2 ^ (A = +,-). (29) 

The index i and j have been omitted in the definition of g^ of Eq. (28) to avoid large chains of 
indices in Section II. The index i always refers to the squark index of flavor q = {w, J, s, c, b} 
of the vertices Fi and r2 and the index j always refers to the stop quark index of and 
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The complex phase of the gluino mass M3 is denoted by 0. The couplings are related by 

Si =gl , ^3 =^4 • 

The unitary squark mixing matrices are given as 

-ft:!:) ■ "'=(1:1:) 

For the squark mass and mixing matrices we use the conventions of Ref. [39]. 
B. Neutralino-squark-quark couplings 

For the neutralino-squark-quark couplings, we again follow the notation of Ref. [39] , where 
explicit expressions for these couplings can be found. For completeness, since the uu- 
channel is the dominant tt production process, we provide here the neutralino-up-quark- 
squark coupling, which reads with the restriction /m,, = GeV [39]: 



5/7 



iVleg^ with (31) 



S4 



S4 



i^-Ki-uf--^-NU-uf_ 

Isw ■ ' 6cvi/ ' ' 

Amws^swNks ■ - 3cwmtNk,4 ■ U]- 



^- {3cwmwspNl2 + swmwspNlij ■ Uj__ + 3cwmtNl,^ ■ Uj _ 



(32) 



where we used the shorthand notations cw = cos 6\y,sw = sin 6w and = sin /3 with tan /3 
— the ratio of 
are related by 



^ the ratio of the two Higgs field vacuum expectation values. Here the coupling parameters 



Si =sT , S3 =sj ■ 

The neutralino mass matrices in the used convention have been taken from Ref. [63]. 



C. Analytic Expressions for the one-loop functions D" 

The partonic differential cross section for the direct box diagrams of Fig. 1(a) and 

crossed box of Fig. 1(b) are given in Eq. (8) and Eq. (10), respectively, in terms of coupling 
parameters and loop functions D". Neglecting the initial-state quark masses and using the 
mass assignments of Fig. 1, the functions D"{s, cos 9) are given with t = mj — j{l— /3f cos 0) 
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as 



++(f,cos0) = ~ 32)^2'- " {—^itntf'Di + (^—mf — mjs + 2mji — t^) mtm2D2] 
D-+++(f,cos0) = -^;^^-mtm4[{-mf-mjs + 2mjf-f)D2-fD3] 

D~++ (f , cos ) = - --^ ■ {Doo ■ (-2mf + 2mjs + 4mjt - if - 4st - if) 

+ rn^fDn + [mfs-mjf - 2mjst + f + 2ft + sf) D13 



+ (mf + m^s- 2mfi+ mj ■ D22 + mjfD23 } 
D"+"+(f,cos0) = -^;^^-m2m4- {m'} + mjs-2mjt + f)Do, (33) 



where = Dijj (O, to^ , to^ , 0, f , 5, m j , , , m|) are written in the convention of Ref . [43] . 
From these expressions the contribution to Fig. 1(b) can be obtained by replacing 
D"(f,cos 6) — )> D"{s, — cos 0) and multiplying by a factor of ( — 1) for exchanging the final- 
state fermions. 
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